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ABSTRACT 

Recently discovered large-scale structure of Lya Emitters (LAEs) raises a novel chal- 
lenge to the cold dark matter (CDM) cosmology. The structure is extended over more 
than 50 Mpc at redshift z — 3.1, and exhibits a considerably weak angular correlation. 
Such properties of LAE distributions appear to be incompatible with the standard bi- 
ased galaxy formation scenario in the CDM cosmology. In this paper, by considering 
the possibility that LAEs are short-lived events, we attempt to build up the picture 
of LAEs concordant with the CDM cosmology. We find that if the lifetime of LAEs 
is as short as (6.7 ± 0.6) x 10 7 yr, the distributions of simulated galaxies successfully 
match the extension and morphology of large-scale structure of LAEs at z = 3.1, and 
also the weak angular correlation function. This result implies that LAEs at z = 3.1 
do not necessarily reside in high density peaks, but tends to be located in less dense 
regions, in a different way from the expectation by the standard biased galaxy forma- 
tion scenario. In addition, we make a prediction for the angular correlation function 
of LAEs at redshifts higher than 3. It is found that the prediction deviates from that 
by the standard biased galaxy formation scenario even at redshifts 4 < z < 6. 
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Recently, the deep imaging surveys by 8 ~ 10 m class tele- 
scopes with narrow-band filter have effectively revealed the 
properties of Lya emit ters (LAEs), which are one class 
of high red shift objects (ICowie fc Hull 19981 ; iHu et alJll99St 
1 1999] , 12002ft . Based on the observational results, it is inferred 
that LAEs have smaller sizes, much less dust, and a smaller 
amount of stellar component than the other class of high red- 
shift gala xies, e.g., Lyman Break Galaxies (LBGs) at same 
redshifts (|Shaplev et all l200ll ; IVenemans etail 120051 ). The 
spatial distribution of observed LAEs generall y shows large 
filamentary structure llShimasaku et al. 1 120031 ; lOuchi et al.l 



12004 120051; iMatsuda et al.ll2005l) . However, as pointed out 



by Hamana et al.l ( 20041 ) recently, the observed properties of 
LAEs such as weak angular correlation function (ACF) are 
not explained well by a standard biased galaxy formation 
scenario in the context of A cold dark matter (ACDM) cos- 
mo logy. Especially, the larg e-scale structure of LAEs found 
by lHavashino et al] l|2004 ) is difficult to reproduce. The 
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large-scale structure shows belt-like structure rather than 
filamentary structure, and may correspond to 6cr density 
fluct uation if it follows und erlying dark matter distribu- 
tion (Kauffma nn et aI1|l9 99). ACF is significantly weaker 
than that predict ed in a conventional biased galax y form a- 
tion model (e.g. iKauffhiann. Nusser fc SteinmetJ l|l997ft ). 
Moreover, ACF becomes negative at small scale of < 
180 arcsec (< 6h~ x M pc) in high density regions (HDR) 
l|Havashino et al.|[20o3 ). Hence, the observational features 
of spatial distribution of LAEs appear to be incompatible 
with the standard biased galaxy formation model. 

From a theoretical point of view, it is recently ar- 
gued that LAEs corresponds to an earl y chemodynamical 
evolut ion phase of primordial galaxies (|Mori fc Umemural 
l2006al lbl). In an ultra high resolution simulation on 
the dynamical and che mical evolution of galaxy by 
iMori fc Umemural i|2006al lbl). it is shown that multiple su- 
pernova explosions at an early phase of < 3 x 10 s yr result 
in forming high density cooling shells, which emit so strong 
Lya as to account for the luminosity of LAEs. However, it 
has not been argued whether this picture of LAEs is consis- 
tent with the observation. 

In this letter, the spatial distributions of LAEs are sim- 
ulated by taking into account the lifetime of the emitters, 
which has not been thitherto considered in the standard 
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Figure 1. Schematic view of the present model. The amplitude 
of density fluctuations at a redshift are shown against the spatial 
location. The density peaks above <5 max , e.g., peak A, have already 
finished their lifetime as LAEs. In contrast, the density peaks 
below <5 m i n e.g., peak C, do not yet start to shine as LAEs. The 
density peaks between <5 m i n and 5 max , e.g., peaks B or D, can 
be observed as LAEs. In this view, peak A is the oldest galaxy 
after LAE phase. Peak B is the youngest LAEs and peak D is the 
oldest LAEs. Peak E has just finished to shine as LAEs at the 
redshift. 

biased galaxy form ation scenario (|KaufImann et al . 1999; 
lHamana et"alll2004 l. Then, we investigate whether the pic- 
ture of short-lived LAE s can explain the clust ering prop- 
erties of LAEs found by lHavashino et alj (|2004j) . In §2, we 
describe the basic picture and numerical method. In §3, the 
results are presented with some discussion. §4 is devoted to 
the summary. Throughout this letter, we adopt ACDM cos- 
mology with the matter density = 0.3, the cosmological 
constant J1a = 0.7, the Hubble constant h = 0.7 in units of 
H = 100 km s " 1 Mpc" 1 , the bary on density Q B h 2 = 0.02, 
and <t 8 = 0.92 jSpergel et al.ll2003n . 



2 MODEL 

2.1 Basic Picture 

In Fig. [TJ the schematic picture of the present galaxy for- 
mation model is presented. In the context of a conventional 
biased galaxy formation model, density peaks with the am- 
plitude which exceeds a minimum threshold value (<5 m i n in 
Fig.[T} in the linear regime are identified as galaxies. In other 
words, only this threshold of fluctuations has been discussed 
as a parameter of biased ga laxy formation (|Kauffmann et ail 
1 1999 ; lHamana et al.|[2004] ). 

Here, we introduce an additional criterion with postu- 
lating that LAEs evolve to galaxies with no strong Lya emis- 
sion after their short lifetime. More specifically, we take fol- 
lowing assumption for LAEs: (i) LAEs are galactic objects 
that form at peaks of density fluctuations, (ii) LAEs are in 
the phase of their first starbursts. (iii) Chemical evolution 
of LAEs results in strong attenuation of Lya emission due 
to the increase of dust, and therefore cannot be observed 
as LAEs after their lifetime. We incorporate this picture by 
setting a maximum threshold of density fluctuations (5 max 



in Fig. [TJ. Then, we regard the fluctuations between (5 m in 
and 5 max as LAEs (a shaded region in Fig. [1} ■ The growth 
time from 5 m in to <5 ma x corresponds to the lifetime of LAEs. 
For instance, peak A in Fig. [1] is the evolved galaxy that 
cannot be observed as LAEs because of exceeding 5 max at 
the redshift. Peaks B and D can be observed as LAEs. Peak 
D is the oldest LAE. 

2.2 Numerical Method 

To compare o ur model with the obs erved clustering proper- 
ties of LAEs l|Havashino et al1l2004h . we numerically gener- 
ate LAE distributions, and estimate the two-point angular 
correlation function by following procedures. 

2.2.1 Generation of LAE Spatial Distribution 

It is assumed that the dynamical evolution of baryonic mat- 
ter follows that of dark matter. Density fields of dark matter 
are created by generating random Gaussian density fields, 
and the dynamical evolut ion is represented by tru ncated 
Zel'dovich approximation (ISathvaprakash et al.lll995l ). This 
approximation traces the growth of density fluctuations in 
the linear regime, and truncate nonlinear growth by sup- 
pressing the amplitude of density fluctuations that becomes 
nonlinear. In the present simulation, we use fc-space Gaus- 
sian window n = exp(— k 2 /2ko) as truncation, where &g 
corresponds to the scale that just enters nonlinear stage at 
a redshift z. The truncated power spectrum of density fluc- 
tuation at z, P*(k,z), is written as 

P*(k,z) = P(k,z)Il 2 (k,z), (1) 

where P(k, z) is the power spectrum of density fluctuations 
at z. The wavenumber ki and real scale n have the relation 

Of fci = 27T /n. 

According to the ACDM theory, the physical size of la 
density fluctuation that collapses just at z — 3.1 is about 
R = 1ft -1 Mpc. In this study, we consider density fluctua- 
tions down to this physical size. In order to directly com- 
pare our model wi th the LAE data in the comoving volume 
of (50/i _1 Mpc) 3 (lHavashino et al.l 120041 ). we simulate the 
same comoving volume with 200 3 grids. The whole simu- 
lation box contains 4.5 x 10 15 Mq in dark matter compo- 
nent, and each cell is (0.25h _1 Mpc) 3 and has 5.7 x 1O 8 M 
on average. Next, we make coarse-graining of density fields 
by comoving volume of {lh" 1 Mpc) 3 which corresponds 
to the physical size of interest. Each coarse-grained cell 
has 3.6 x 10 10 Mq in dark matter on average. The coarse- 
grained cells that satisfy the density fluctuation criterion, 
5mm ^ 5 ^ <5 m ax, are regarded as LAEs. The positions are 
assumed to be the center of mass in coarse-grained cell. 

We choose several combinations of <5 m i n and <5 max . A set 
of 5 m i n and <5 max is constrained so that the number of sim- 
ulated LAEs should match the obs erved number of LAEs 
at z = 3.1 (|Havashino et all 120041 ) . Thus, if <5 m m is set, 
then <5 max is determined from the constraint of the num- 
ber of LAEs. In the linear regime, only density flu ctuations 
with 8 ^ 1.7 corresponds to collapsed objects (|Peacockl 
1998). Therefore, we consider <5 m j n larger than 1.7. Resul- 
tant three-dimensional distributions of LAEs are projected 
into a two-dimensional plane to compare observed angular 
distributions. 
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Figure 2. LAE distribution for different values of <5 m i n and <5 ma x- Left panel is the LAE distribution for the model with <5 m i n = 1.7 and 
5max = 1.75, and middle panel is <5 m j n = 2.5 and <5 max = 2.63. Right panel is <5 m j n = 5.4 and <5 max = oo, which is corresponding to a 
conventional bia s ed ga laxy formation model. The contours represent high density regions (HDRs) of LAEs under the same condition as 
Havashin o et al.1 II2004T) 



2.2.2 Angular Correlation Function (ACF) 

To calculate the two-point ACF of the simulated spatial dis- 
tribution of LAEs, we use a following well-known estimator, 



w(9) 



N T (DD(8)) 



(2) 



N s (DR(9)) ' 

l|Peeblesll980l ; |Peacocklll99Sn , where N s and N v are the mean 
surface number density of simulated LAEs and that of ran- 
domly distributed points (RDPs), respectively. RDPs are 
distributed over the same area as LAEs. (DD{9)) is the av- 
eraged pair number of LAEs in a range of (9, 9 + d6), and 
(DR(9)) is the averaged pair number between LAEs and 
RDPs in a range of (9,9 + d9). To raise the precision of 
statistics, we calculate ACFs for 30 different realizations of 
density fluctuations and average them. Then, the error on 
w(9) is defined by the standard deviation of ACFs. 



3 RESULTS 

3.1 Clustering Properties of LAEs 

In Fig. [5J the spatial distributions of simulated LAEs are 
shown for different values of <5 m i n . The contours depict 
'High D ensity Region' (HDR) d efined under the same con- 
dition as lHavashino et al.l (|2004l ) , where the number density 
smoothed with a Gaussian kernel of ctg = 90 arcsec (corre- 
sponding to 3/i -1 Mpc) is equal to the mean number density 
in the entire field. The left panel in Fig. [2] is the model with 
<5min = 1-7 and (S m ax = 1.75, and the middle panel is the 
model with <5 m i n = 2.5 and 5 max = 2.63. In the right panel, a 
conventional biased galaxy formation model is shown, where 
Smln ~ 5.4 is assumed and all fluctuations with 8 ^ <5 m i n 
are regarded as LAEs. In all these panels, we can recognize 
large-scale structures, but the clustering manners are some- 
what different. The spatial distributions in the biased galaxy 
formation model exhibit very strong contrast and their clus- 
tered regions are fairly isolated. On the other hand, the dis- 
tributions for (S m iu, <5max) = (1.7, 1.75) or (2.5, 2.63) appear 



to be belt-like and less clu stered, similar to the ob served spa- 
tial distribution of LAEs l|Havashino et al.|[200i) . 

In order to quantify the difference in spatial distribu- 
tions, we calculate ACFs. In Fig. [3] the resultant ACFs for 
all the model are pres ented. Also, the ACF for LAEs ob- 
served in SSA22a field l|Havashino et al.|[2003 ') is shown. The 
upper and the lower panels show the ACF in the whole 
region and that in HDR, respectively. The results show 
different behaviors on scales smaller than ~ 300 arcsec. 
The biased galaxy formation model shows strong corre- 
lation on small scales as expe cted in a standard biased 
model (jKauffmann et all 1 19991 ). and obviously does not 
match the ACF of observed LAEs in SSA22a field. Fur- 
thermore, the model with (<5 m i n , <5max) = (2.5, 2.63) results 
in slightly stronger ACF than the observation. The model 
with (<5 m in, <5 max ) = (1.7, 1.75) remarkably agrees with the 
ACF of LAEs in SSA22a field. In the HDR, the ACF exhibits 
negative correlation in the same way as the observation. The 
reduction of ACF for smaller (5 m i n is understood as follows. 
In the random Gaussian density fields in a ACDM universe, 
higher density peaks are more clustered, while lower density 
peaks are located in less dense regions surrounding highest 
density regions. Thus, if small <5 m j n is adopted and high- 
est peaks are cut by <5 ma x, the objects of interest are lo- 
cated in less dense regions and accordingly the amplitude of 
ACF becomes smaller. Hence, the result that the model with 
(<5min, <5max) = (1-7, 1.75) reproduces the observed ACF im- 
plies that LAEs at z = 3.1 do not reside in highest density 
peaks, but are located in less dense regions. 

Observationally, LAEs have been discovered around 
known overdensities that generally indicate strong correla- 
tion such as proto-clus ter region including massive galaxies 
such as radio galaxies ( Steidel et al J 12000 : lHavashino et al.1 
|2004 IVenemans et al.1 |2005T )~ The observed overdens ties 
may correspond to the situation shown in the right pane 
of Fig. [2j and the observed LAEs correspond to the left 
panel of Fig. [5] In that sense, the results here look consis- 
tent with these observational features. Hence, the picture 
in this paper can explain not only a correlation function 
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Figure 3. Two-point angular correlation function (ACF) of LAE 
distributions for each model. Upper and lower panels show ACFs 
in the whole region and in HDR, respectively. Open squares rep- 
resent ACFs for 

^min — 1-7 and 5 m ax — 1.75, and open circles 
for 5 m ; n = 2.5 and 5 max = 2.63. Open triangles represent ACFs 
for ^min — 5.4 and <5 max — oo, which is corresponding to a biased 
galaxy f ormation model. Cross es are ACFs of LAEs observed in 
SSA22a (Havas hino et al.ll2004n . Note that the upper panel is dif- 
ferent from the lower panel in the scale of vertical axis. 



but also other clustering properties of LAEs such as the 
morphology of HDR and the environment where LAEs at 
z ~ 3 is discovered. According to the recent large survey 
such as SDSS, ACFs of late-type galaxies show weaker cor- 
relation compared w ith that of early- type galaxies at z < 0.1 
|Zehavi et all [2002). That is to say, as well known, late- 
type galaxies are located at lower density fields. As shown 
here, LAEs at z = 3.1 should be located in less dense re- 
gions. Hence, it is suggested that a large fraction of LAEs 
at z — 3.1 may be the precursors of late-type galaxies. 



3.2 Lifetime of LAEs 

As shown above, 8 m in = 1-7 gives the best fit model to ac- 
count the observed ACF. Si nce 5 = 1.7 is a critical amplitude 
for a fluctuation to collapse l|Peacocklll998h . we can conclude 
that LAEs begin to shine just after the collapse. In other 
words, LAEs should be in the first phase of galaxy evolution. 
Since the model with (<5 m i n , 5 max ) = (1.7, 1.75) agrees with 
the observed ACF at z = 3.1, LAEs are thought to shine 
during the growth time from 5 m in to 5 max . The fluctuation 
with <5 max at z = 3.1 collapses at a higher redshift z co \\ when 
the amplitude exceeds <5 m i n . Hence, the lifetime of LAEs can 
be assessed by the cosmic time between z = 3.1 and z co ii, 
which is 6.7 x 10 7 yr. Here, there is a small uncertainty in 
this estimation. When <5 m i n is chosen, <5 max is determined 
to match the number of observed LAEs. Since we generate 
random numbers to produce density fluctuations, a different 
set of random numbers results in slight difference in <5 max . 
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Figure 4. Two-point angular correlation function (ACF) of sim- 
ulated LAEs at redshifts at z = 3.1, 4.0, 5.0, and 6.0. Open circles 
are short-lived LAE model, while filled circles are biased galaxy 
formation models. In the short-lived LAE model, density peaks 
between <5 m i n = 1.7 and <5 m ax = 1.75 are regarded as LAEs. In the 
biased galaxy formation model, the number of objects is scaled 
as to be the same as that in the short-lived LAE model. So, 5 m i n 
is set to 5.4, 4.2, 3.5, and 2.8, at z = 3.1, z = 4, z = 5, and z = 6, 
respectively. 



For the model of <5 m i n = 1.7, we have <5 max = 1-75 ± 0.01 
as a result of 30 different realizations. Then, the lifetime of 
LAEs is estimated to be (6.7 ± 0.6) x 10 7 yr. Similarly, for 



the model of S n 



2.5, we have <5 max = 2.63 ± 0.02. Then, 



the lifetime is slightly longer as (2.0 ± 0.4) x 10 yr. 

This result on LAE lifetime nicely agrees with an upper 
limit that is argu ed by realistic numerical simulations for 
galactic evolution l|Mori fc Umemu ra 2006a||3). 



3.3 Luminosities of LAEs 

We also calculate Lya luminosities of simulated 
LAEs, using an evolutionary spectra l synth esis code 
'PEGASE' l|Fioc fc Rocca-Volmerange] Il9971 ). As a 
result, we have found that evaluated Lya lumi- 
nosities match those of observed LAEs (Ll V c» ~ 



10 42 ~ 4 3 ergs~ 1 ) dHavashino et all 12004 iMatsuda et all 



|2004| ; Ivan Breukelen. Jarvis. fc Venemand l2005h . In this 
paper, density fields are coarse-grained by a scale of 
l/i -1 Mpc which corresponds to la density fluctuations 
in the ACDM cosmology. If a smaller scale is taken, 
intrinsic Lya luminosities fall short of 10 42 ergs _1 during 
Lya bright phase. For instance, if a coarse- graining scale is 
0.25/i -1 Mpc, intrinsic Lya luminosities are ~ 10 41 ergs -1 . 
On the other hand, if a scale larger than lfr -1 Mpc is taken, 
the number of collapsed objects is not enough to accout for 
the observed LAE number. Hence, lcr density fluctuations 
are favorable to explain the observations. 



3.4 ACF of LAEs at 3 < z < 



By assuming the best fit model (<5 m i n = 1.7 and <5 ma x 
1.75), we can predict ACFs of LAEs at higher redshifts. 



In 
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Fig. H the prediction of ACFs at redshifts of 3.1, 4.0, 5.0, 
and 6.0 are presented. A biased galaxy formation model is 
also presented, where the number of objects is scaled as to be 
the same as that in the best fit model. 5 m in of biased galaxy 
formation model in each redshift is <5 m in = 4.2 at z = 4, 
5 m in = 3.5 at z = 5, and <5 m i n = 2.8 at z = 6, respectively. 
As seen in this figure, the ACF of best fit model approaches 
to that of biased galaxy formation model at higher redshifts. 
In other words, a larger fraction of collapsed objects shine as 
LAEs at higher redshifts. But, it is worth noting that there 
is still noticeable difference between the best fit model and 
a biased galaxy formation model even at z = 6. It implies 
that a certain fraction has been already extinguished, so that 
they are not detected as LAEs. 



4 SUMMARY 

To account for the rec ently discovered large-scale structure 
of LAEs at z = 3.1 (|Havashino et al1l2004T ). we have in- 
troduced a novel picture for LAEs by focusing on the life- 
time of emitters. We have simulated the spatial distribu- 
tions of collapsed objects by generating random Gaussian 
fluctuations based on the truncated Zel'dovich approxima- 
tion in the ACDM cosmology. We have found that a conven- 
tional biased galaxy formation model is not reconciled with 
the observed correlation function of LAEs. If highest peaks 
above S = 1.75 are cut and mild peaks between S = 1.7 
and 5 = 1.75 are regarded as LAEs, the clustering proper- 
ties including two-point angular correlation function agree 
quite well with the observation. Lya luminosities also match 
those of observed LAEs. The growth time from 8 — 1.7 to 
8 — 1.75 can be translated into the lifetime of LAEs, which 
is assessed to be (6.7 ± 0.6) x 10 7 yr. A fluctuation with 
8 — 1.7 corresponds to an object that just collapses at the 
redshift. Thus, LAEs are thought to be in the early evo- 
lutionary phase o f galaxies, consistently wi th a recent theo- 
retical prediction dMori fc Umemur a 2006a b]). We have also 
predicted the correlation function at redshift higher than 3 
in the picture of short-lived LAEs. It is suggested that a cer- 
tain fraction of young galaxies have already ended the LAE 
phase even at redshift z — 6. 
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